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ABSTRACT 

Scattering of Lya photons by neutral hydrogen gas in a single outflowing 'supershell' 
around star forming regions often explains the shape and offset of the observed Lya 
emission line from galaxies. We compute the radiation pressure that is exerted by this 
scattered Lya radiation on the outflowing material. We show that for reasonable phys- 
ical parameters, Lya radiation pressure alone can accelerate supershells to velocities 
in the range Ush = 200-400 km s~^. These supershells possibly escape from the grav- 
itational potential well of their host galaxies and contribute to the enrichment of the 
intergalactic medium. We compute the physical properties of expanding supershells 
that are likely to be present in a sample of known high-redshift (z — 2.7-5.0) galaxies, 
under the assumption that they are driven predominantly by Lya radiation pressure. 
We predict ranges of radii Tgh = 0.1-10 kpc, ages tsh = 1-100 Myr, and energies 
Esh = lO^'^-lO^^ ergs, which are in reasonable agreement with the properties of local 
galactic supershells. Furthermore, we find that the radius, rgh, of a Lya-driven super- 
shell of constant mass depends uniquely on the intrinsic Lya luminosity of the galaxy, 
Lq, the HI column density of the supershell, A^hi, and the shell speed, Wsh, through 
the scaling relation Tgh oc La/ {N-mvl-^- We derive mass outflow rates in supershells 
that reach ^ 10-100% of the star formation rates of their host galaxies. 

Key words: cosmology: theory-galaxies: high redshift-radiation mechanisms: 
general-radiative transfer-ISM: bubbles 
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1 INTRODUCTION 

The Lya emission line of galaxies is often redshifted rel- 
ative to metal abs orption lines (e.g. iPettini et al.l l200ll : 
IShaplev et al.l [20031 ). Scattering of Lya photons by neutral 
hydrogen atoms in an outflowing 'supershell' surrounding 
the star forming regions can naturally explain this obser- 
vation, as well as the typically obser ved asymmetric spec- 
tral shape o f the Lya emission line llLequeux et al.l 19951: 
Lee fc AhnI Il 998: Te norio-Tagle et all Il999l: lAhn fc Led 



atomic hydrogen around HII regions i s confirmed by HI- 
observations of our own (iHeiled [l984l ) and other nearby 
galaxies (e.g. iRvder et al.lll995l ). The largest of these shells. 



so-called 'supershells', have radii of r 



1 kpc (e.g 



iRvder et al.Hl995l : iMcClure-Griffiths et al.ll2002l . |2006| ) and 

HI column densities in th e range A^hi ~ lO^ ^ -lO^'^ cm 



e.g. L cqucux ot all 1 19951 : iKunth et ahl Il998l : lAhnl |2004| : 



20021: 



Ahn et al 



Verhamme et al. 



.2003i: 



l2006l , l2008l ). The presence of an outflow 

may also be required to avoid complete destructiorQ of the 
Lya radi ation by dust and to allow its escape f rom the host 
galaxies llKunth et al.l l l998l : iHaves et al]|2008l : lOsthn et af] 



I Verhamme et al.ll2008D . The existence of larger extragalac- 
tic supershells has been inferred from spatially extended 

m" ^jj^^^ '^'^t"^ l|'' 2003l 1 AhnI bocS' ^ ^'P'^) Lya P-cygni profi les, that were observe d 

^ — esse_e — ^, J, around local star burst galaxies l|Mas-Hesse et all |2003| ). 

Thes e 'supershells' differ f rom more en ergetic 'superwinds' 
(e.g. iHeckman et aT]|l990l : lMartinll2005l ). in which galactic- 
scale biconical outflows break-out of the galaxies' interstel- 
lar medium with high velocities (up to ~ 10^ km s~^, of 
which M82 is a classical example). Supershells are thought 
to be generated by stellar winds or supernovae explosions 
whi ch sweep-up gas into a thin expandin g neutral shell (see 
e.g. iTenorio-Tagle fc Bodenheimerlll988l . for a review). The 
back-scattering mechanism attributes both the redshift and 
asymmetry of the Lya line to the Doppler boost that Lya 
photons undergo as they scatter off the outflow on the far 
side of the galaxy back towards the observer. 

In this paper we show that the radiation pressure ex- 



l2008l : lAtek et al.ll200a ). 

Indeed, the existence of outflowing thin shells (with 
a thickness much smaller than their radius) of neutral 



* E-mail:mdijkstr@cfa. harvard.edu 

^ Lya may also avoid complete destruction by interstellar dust 
when it is primarily locked up in col d clu mps embedde d in a hot 
medium that is transparent to Lya l|Neufeld,1991. : .Hansen fc OhI 
|2006|) . 
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erted by backscattered Lya radiation on the outflowing gas 
shell can accelerate it to velocities that reach hundreds of km 
s~^. Our basic result can be illustrated through an order- 
of-magnitude estimate. The net outward force exerted on 
a shell of ga^ by backscattered Lya radiation of luminos- 
ity Lbs, a is d(msh«sh)/rfi = -I/BS.a/c. Here, rrish is the to- 
tal mass in the shell, Vsh is the shell velocity and c is the 
speed of light. If we assume that backscattering occurs over 
a timescale tss, and the shell mass is constant in time, 
then the net velocity gain by the shell is Avsh ~ 250 km 
s-i (m,h/10^ Mo)-'i(fBs/50 Myr)(LBS,c./10*3 erg/s). The 
adopted rrish value corresponds to a thin spherical HI shell 
with a column density A^hi = 10^" cm~^ and radius r = 1 
kpc. The fiducial values of Lbs, a and tss were chosen to 
yield the typical observed shell speed and exp ected shell life- 
time in Lyman-brea k galaxies (LBGs; see e.g. lShaplev et all 
l2003l : IVerhamme et al. 2008). Here, and throughout this pa- 
per, we assume that the supershells are mostly neutral. This 
assumption is based on 21-cm observations of local super- 
shells which indicate that these are very thin (thickness <^ 
their radius, e.g Heiles 1984) and therefore dense. At these 
high densities one expects free electrons and protons to re- 
combine efficiently and make the gas significantly neutral. 
However, iPettini et all ||2000D argue that, based on the ob- 
served reddening of the UV continuum spectrum, either the 
dust-to-gas ratio in the supershell of LBG MS1512-cB58 is 
higher than the galactic value, or alternatively that the HI 
only makes up /hi ~ 14 — 33% of the total gas mass of the 
shell, with the remainder of the gas being molecular or ion- 
ized. We discuss how this latter possibility affects our results 

in m 

This paper is a closely related to another paper in which 
detailed Monte-Carlo Lya radiative transfer calculations are 
performed to assess Lya radiation pressure in a more general 
context (Dijkstra & Loeb 2008, hereafter Paper I). These 
radiative transfer calculations illustrate clearly that Lya ra- 
diation pressure may be important in driving outflows in 
various environments. Importantly, the pressure exerted by 
Lya radiation alone can exceed the maximum possible pres- 
sure exerted by continuum radiation through dust opacity 
(see Paper I), which in turn can exc eed the total kineti c 
pressure exerted by supernova ejecta (|Murrav et al.|[2005h . 
This implies that Lya radiation pressure may thus in some 
cases provide the dominant source of pressure on neutral 
hydrogen gas in the interstellar medium (ISM). 

The goals of this paper are (i) to gauge the general im- 
portance of Lya radiation in the particular context of out- 
flowing galactic supershells, and (ii) to explore the related 
observable properties of Lya-driven supershells. In § [5] we 
compute the time evolution of the velocity of an expanding 

^ The net momentum transfer rate per photon that enters the 
shell is Ap = — ii)er, where denotes a unit vector that 

is pointing radially outward. Furthermore, fi = kjn • kout, where 
kin (kout ) denotes the propagation direction of the Lya photon as 
it enters (leaves) the shell. When averaged over all directions, the 
total momentum transfer rate is Na^^^^Sr, where Na is the rate 
at which Lya photons enter the shell. This total net momentum 
transfer rate applies regardless of the number of times each Lya 
photon scatters inside the shell, or regardless of whether the Lya 
photon is absorbed by a dust grain (and possibly re-emitted in 
the infrared). 



shell of neutral gas that encloses a central Lya source. We 
consider a range of models, which are subsequently applied 
to known Lya emitting galaxies at high redshifts. Finally, 
we discuss our results and their implications in § [3] The 
parameters for the background cosmology used throughout 
our discussion are (n^, Oa, O;,, fe) = (0-27 0.73, 0.042, 0.70) 
(|Dunklev et al.ll2008l : Komatsu et al.ll2008l ). 



2 PRESSURE FROM BACKSCATTERED LYq 
RADIATION 

As mentioned in § [1] the net momentum transfer rate from 
the Lya radiation held to the shell is given by 

dt ^ '"^^ = /--(-^l^'^Hl) — , (1) 

where La is the total Lya luminosity of the galaxy, and 
/scat(^'sh, A^Hi) denotes the fraction of Lya photons that are 
scattered in the shell. The latter depends on various factors, 
including the column density of neutral hydrogen atoms in 
the shell, the intrinsic Lya spectrum (i.e. the emitted spec- 
trum prior to scattering in the shell), as well as the shell 
velocity. For example, as Uah increases, atoms in the shell in- 
teract with Lya photons that are farther in the wing of the 
line profile, and /scat diminishes. Eq. [T] assumes that the HI 
shell mass remains constant in time. This assumption is mo- 
tivated by the observation that supershells are well defined 
shells of HI gas that are physically separate from the central 
Lya source. These outflows differ from steady-state outflows 
(i.e. outflows in which mass is ejected at a constant rate 
M) that occur aroun d late-type stars fe.g. ISalpeteij Il973 : 
llvezic fc Eiit^liggsl ). In reality, an HI supershell is not 
likely to maintain a constant mass as it sweeps through the 
ISM, but its evolution depends on the detailed hydrodynam- 
ics of the shell as it propagates through the ISM, which in 
turn depends on the assumed properties of the ISM. In this 
first study, we focus on simplified models that gauge the gen- 
eral importance of Lya radiation pressure on the dynamics 
of galactic supershells. 

Furthermore, our calculation of the parameter /scat con- 
servatively ignores that possibility of 'trapped' Lya photons 
by the expanding supershell, which can boost the radiation 
pressure considerably. If radiation trapping were included, 
then the parameter /scat would equal the 'force-multiplier' 
Mf, which was shown to be much greater than unity at low 
Wsh in Paper I (see § [3] for a more detailed discussion) . We 
adopt a conservative approach because the calculation of 
Mf in Paper I assumed that the expanding shell contained 
no dust. Under realistic circumstances, the existence of dust 
inside the supershell will suppress the number of times the 
Lya photons can 'bounce' back and forth between opposite 
sites of the shell, which leads to a reduction in the value 
of Mf- However, provided that the dust is located inside 
(or exterior to) the HI shell, the true value of Mp is always 
larger than the /scat adopted in the paper (see §[3|. 

2.1 Calculation of /scat 

We compute /scat(i'sh, A^'hi) under the conventional set of 
assumptions: (i) the emitted Lya spectrum prior to scat- 
tering is assumed to be a Gaussian with a velocity width 
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Table 1. Parameters of Models used in Figure [2] 
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Figure 1. The fraction of scattered Lya photons /scat as a func- 
tion of shell velocity Vc^i^ for a range of column densities Afjji (in 
cm~-^, for (T = 75 km s~^) and emitted line widths quantified by 
a (labeled by 75, 250 and 500 km s"! for logAfni = 20). Higher 
A^Hi values require larger shell speeds in order for the Lya pho- 
tons to propagate through the shell unobscured, because as A^hi 
increases the photons need to be farther in the wing of the line 
profile (in the frame of the shell) for them not to scatter. Fur- 
thermore, a larger intrinsic line width, c, serves to flatten the 
dependence of /scat on the shell speed (see text). 



a = Ucirc (jSantosll200l : iDiikstra et al.ll2007l ). where Wdrc is 
the circular virial velocity of the host dark matter halo; and 
(ii) the outflow is modeled as a single expanding thin shell of 
gas with a column density A^hi (a s in lAhn et aDl2003l : lAhnl 
l2004l : IVerhamme et al.ll2006l . [ 20081 ) . We can then write 



10^° cm 
scat increases with increasing 



/scat(wsh, A'"hi) = ^= — / dx e 

(2) 

where frequency is denoted by the dimensionless variable 
X = {v — i'o)/Auo, with Avn — Vthi^o/c and Vth being the 
thermal velocity of the hydrogen atoms in the gas given by 
Vth ~ \/2kBT/mp, ks is the Boltzmann constant, T — lO** 
K is the gas temperature, irip is the proton mass, i/g = 2.47 x 
10^^ Hz is the Lya resonance frequency, ax = cr/iith, and 
Xah = I'sh/wth. In Figure[T]we show /scat as a function of Ush 
for a range of column densities A''hi (for a = 75 km s~^) and 
emitted line widths a (for A^'hi 

As illustrated in Figure [T] /. 
A'hi, because photons need to be farther in the wing of the 
line profile (in the frame of the shell) for them not to scatter. 
For example, Figure[l]shows that ~ 50% of the Lya photons 
are scattered when the shell speed is ~ 200 km s"^ for A*'hi = 
10^9 cm~^, whereas the same fraction is reached when the 
sheU speed is ~ 600 (2000) km s'^ for TVhi = 10^" (iVni = 
lO'^^) cm"'^. Furthermore, increasing the intrinsic line width 
serves to fiatten the dependence of /scat on the shell speed, 
because at large a and low shell speeds there is a substantial 
fraction of photons far in the wing of the line profile. 



# 




La 


''sh.i 


fsh.i 






(km s-l) 


(erg s-i) 


(kpc) 


(km s-i) 


(cm~^) 


1. 


150 


10« 


1.0 


50 


1020 


2. 


150 


10« 


1.0 


150 


1020 


3. 


150 


10« 


1.0 


50 


lOlo 



2.2 The Time Evolution of the Shell Speed: 
General Calculations 

Given /scat(wsh, A^hi) we next use equation H]) to compute 
the evolution of the shell position (rsh) and speed (Vsh) as 
follows; 

1. The initial location and shell speed are denoted by rsh,i 
and Vsh,i- If the outflowing shell covers 4-k stera dians of the 
sky surrounding the Lya emitting region (see lAhnl |2004| : 
IVerhamme et aL I I2OO8I . for more detailed discussions), and if 
all the shell material is in neutral atomic hydrogen gas then 
msh = 47rrsij [A'hi,! Tip, where r is the radius of the shell, 
nip = 1.6 X 10~^* g is the proton mass, and A^hi.i is the 
initial column density of neutral hydrogen atoms. 

2. We combine equations ^ and to compute dvs\i/dt. 

3. Over an infinitesimal time step At a new shell position 
and velocity are then generated as rsh(t + At) — rsh(t) -I- 
Vsh{t)A.t and Va^{t + At) = i;sh(i) + (dUsh/rft)At. At the 
new position, the shell's column density reduces to A^hi = 
A'Hi,i(T'sh/''sh,i)~^, aud we return to step 2. 

Our fiducial model has the following parameters. We 
adopt a = 150 km s~^, corresponding to the circular virial 
velocity of a dark matter halo of mass M = 10" M© at 
z — 5.7 as appropria te for a typical hos t of known Lya 
emitting galaxies (e.g. iDiikstra et al. I l2007h . We also adopt 
La — 10*^ erg s~i, corresponding to t he typical observe d 
luminosity of Lya emitting galaxies (e.g. lOuchi et al.ll2007l ). 
This Lya luminosity corresponds to a star formation rate 
of ^ 5Mq yr~^ f or a Salpeter I MF and a gas metallic- 
ity Z = O.OSZq (Sch aereil [200l . for Z = Zq this star 
formation rate is higher by a factor of 2). For compar- 
ison, the total kinetic luminosity in supernova ejecta is 
Ln,ech ~ 2xl0*2(£^,,h/10" erg)(A/'sN/10-2) erg s'S where 
A/sN is the supernova rate per unit star formation rate, and 
Emcch is the t otal kinetic energy in ejecta per supernova 
explosion (e.g. iMurrav et al]|2005l . and references therein). 

For the shell, we assume A/hi,! ~ 10^" cm"^^ j-^j^ j — i.Q 
kpc (see !JT] and Verhamme et al. 2008) and Ush,i = 50 km 
s^i. We find that the final shell speed is not very sensitive 
to Wsh.i. The above shell parameters imply a shell mass of 
msh ~ lO'^ Mq. The parameters of the fiducial model are 
summarized in Table [T] 

The time evolution of the shell speed in the fiducial 
model is plotted as the solid line in the left panel of Fig- 
ure [21 which shows that Lya scattering accelerates the shell 
to ~ 100 km s~i after 10 Myr, consistently with the crude 
estimate given in § [T] After 100 Myr, the shell speed reaches 
~ 230 km s~i. The acceleration of the decreases with time 
because /scat decreases with time. As the shell accelerates 
and expands, the fraction of Lya photons that are scattered 
in the shell and the associated momentum transfer rate to 
the shell, decrease. Other lines in this panel represent models 
in which one of the model parameters was changed. 
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Figure 2. Left panel: Time evolution of the speed of an outflowing shell of neutral hydrogen that is accelerated by Lya radiation from 
a star- forming galaxy. The solid line represents our fiducial model (see text and Table [!}, in which hya scattering accelerates the shell 
to ~ 100 km s-l in 10 Myr and to ~ 230 km s'^ in 100 Myr. Other curves represent models in which one of the model parameters was 
modified (except for the specific example illustrated by the blue dotted line; see 12.31 1. The curves demonstrate that scattering of hya 
photons off an expanding shell can accelerate the shell to speeds that reach a few hundred km s~^. Central & Right panels shows the 
shell velocity as a function of its radius & column density, respectively. In all examples except the dotted line, the shell starts at r = 1 
kpc with a modest column density. In the model of LBG MS1512-cB58 (dotted line), the shell starts at r = lO"'^ kpc with a higher 
column density that declines to the observed value as the shell expands (see text). 



The red dashed line shows the shell evolution for a model 
in which the initial shell velocity was increased to Vsh,i = 150 
km s~^. The difference in shell speed between this case and 
the fiducial model decreases with time. In the model repre- 
sented by the green dot-dashed line the shell mass is reduced 
to msh = W^Mq by reducing A^Hi,i. The shell speed evolves 
much faster, and reaches 200 km s~^ after ~ 5 Myr. The 
final shell speed is ~ 400 km s^^. The central (right) panel 
shows the shell velocity as a function of its radius (column 
density). In the fiducial model (and those in which one pa- 
rameter was changed), the shell started out at r = 1 kpc 
with a relatively low column density. The model represented 
by the blue dotted lines is discussed in § 12.31 

In aU the above models, Lya radiation pressure acceler- 
ates the shell up to a few hundreds of km s"'^. For compari- 
son, the escape velocity from a halo with Vdic = 150 km s~^ 
(which corresponds to a mass of 3 X 10"[(l + 2)/4]-=^''2^Q, 
e.g. Barkana & Loeb 2001) for a shell starting at r = 1 kpc 
is VcBc ~ 440 km s"'^ (at 2 = 3, with a very weak redshift 
dependence, see § I2.3|l . Therefore, it is reasonable to claim 
that shells that are driven by Lya radiation pressure may 
escape from massive dark matter halos. Next we consider 
whether this mechanism could operate in known observed 
high-redshift LBGs that have measured values of La, Nm 
and Vsh- 



2.3 The Time Evolution of the Shell Speed in 
Known Galaxies 

In real galaxies, Lya radiation is not the only process that 
determines the shell kinematics, and Equation ((TJ modifies 
to 

rUsh^^^ = /scat(l'sh, Nm)— ~ -Fgravlr) + Attt^ AP, (3) 

dt c 

where f'grav('") is the gravitational force on the shell, and 
AP = Pint — Pcxt with P — pCg being the pressure of 



the medium inside ("int") or outside ("ext") of the shell 
(Elmcgrccn fc Chiang [1983 ). Here, p is the density of the 
medium, and cs is its sound speed. The second and third 
terms on the right hand side of Eq. ([3]) require assumptions 
about the unknown shape of the gravitational potential near 
the star forming region and the local properties of the sur- 
rounding ISM. 

To first order, gravity can be ignored in those galaxies 
for which the calculated terminal shell speed exceeds the 
escape velocity of the dark matter halo. More specifically, 
gravity can be ignored if the Lya momentum transfer rate 
exceeds the gravitational force, La/c ^GM{< rsh)?Tish/'~shi 
where M(< rsh) is the total mass enclosed by the supershell. 
Quantitatively, Lya radiation pressure exceeds gravity when 
the Lya luminosity exceeds 



L.^10«ergs-^^«'^^'^)^^ "^^'^ 



V lO«Af0 / V 106 Mq / V 0.1 kpc 



(4) 



If the radial density profile can be modeled as an isothermal 
sphere, then M(< rsh) ~ Mtot{rsh/rvir) and the required lu- 
minosity increases at smaller shell radii. However, if the shell 
starts at a small radius (i.e. rsh.i ~ 0.01 kpc), then for a fixed 
shell mass the column density is high (e.g. A^hi.i 3> 10^^ 
cm~'^). This, in combination with the small initial veloc- 
ity of the shell, implies that the Lya photons are efficiently 
trapped and the impact of radiation pressure is bigger than 
estimated above (see Fig. [3] and Paper I). In Paper I we 
showed that in this regime, a simple order-of-magnitude es- 
timate for the importance of radiat ion pressure is obtained 
based on energy considerations (e.g. ICox|[l985l : lBithelllll990l : 
lOh fc Hairnanll2002l ). Initially, the total gravitational bind- 



ing energy of the shell is \Ush\ ~ GM{< rah,i)?7iah/rsh,i- The 
energy in the Lya radiation field is Ua = Lattiap, where 
ttrap is the typical trapping time of a Lya photon in the 
neutral shell. The trapping time scales with the shell op- 
tical depth at the Lya line c enter, tq, through the rela- 
tion ttrap ~ 15(To/10^-^)rsh,i/c (lAdamslllOTSl : iBonilha et all 
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1 19791 ). Radiation pressure unbinds the gas when Ua > | Ush \ ■ 
We find this to be the case for all galaxies discussed below 
(except the first model of FDF 1267). Also Uc > \Ush\ for 
those galaxies in which the calculated terminal wind speed 
exceeds the escape velocity of the halo, which confirms that 
ignoring gravity is justified best for those galaxies in which 
the Lya radiation pressure can accelerate the shell to a speed 
that exceeds the escape velocity of the host dark matter halo. 

Since we investigate the properties of Lya-driven super- 
shells, we assume that AP is less than the Lya radiation 
pressure (this assumption is motivated by the possibility 
that Lya radiation pressure provides the dominant source 
of pressure on neutral hydrogen gas in the ISM, see §[T]and 
Paper I) . Hence, we solve Eq. [3] under the assumption that 
the second and third terms on its right hand side can be 
ignored. We initiate a shell at rest {vsh,i = 0) at a very small 
radius (rsh,i ~ 0.01 kpc, and not at the origin to avoid a 
divergent A^Hi,i — > oo). We get constraints on the mass, age 
and radius of the supershell by requiring that it reaches its 
observed expansion velocity at its observed column density 
(see l2.3.1l for an example). 

Our results are summarized in Table (2] The first col- 
umn contains the name of the galaxy. The second column 
contains the HI column density and the third contains the 
shell speed as derived from fitting to the observed Lya 
fine profile ijSchaerer fc Verhammd I2OO8I : IVerhamme et all 
|2008| ). The fourth column contains the intrinsic Lya lumi- 
nosity, which was inferred from the UV-based star forma- 
tion rate, and the intrinsic Lya eq uivalent width der i ved b y 
Schaerer & Verhamme (2008) and I Verhamme et all ||2008D . 
aE ~ (SFR/Mo/yr) x (EWint/70 A) x lO''^ erg s'^ (e.g. 
iDiikstra fc Wvithe 2007 ) . The radius {rsh), angular scale on 
the sky (da), mass (jTigh), and age (tsh) of the supershell are 
given in the fifth, sixth, seventh, and eighth columns, respec- 
tively. The ninth column contains the energy of the super- 
shell, which is given by E^h = ^rnehVah- The tenth column 
gives the asymptotic (terminal) shell speed, to be compared 
with the escape velocity of the galaxy in the eleventh col- 
umn. For an isothermal sphere, the escape velocity relates 
to the circular velocity as Vcbc ~ V^VciTc\/ln[rvir /rsh.i] (see 
Eqs. (2-25), (4-117) & (4-123) of Binney & Tremaine 1987, 
and Eqs. (24) & (25) of Barkana & Loeb 2001). The circu- 
lar velocity Hcirc was obtained from iicirc=FWHM/2.35, in 
which 'FWHM' is the intrinsic Full Width at Half Maxi- 
mum of the Lya line quoted by Verhamme et al. (2008, see 
e.g. Dijkstra et al. 2007 for a more extended discussion on 
this relation). The twelfth column contains a fudge factor / 
which is defined as / = '"sh/(|ashSh)5 where Ssh = La/cmsh, 
and ?ah = 2rsh/i;sh (i.e. / = 1 for /scat = 1, and / < 1 when 
/scat < 1). The radius of the shell can then be expressed as 
a function of the observed quantities La, Nm, and Vsh as 

r-sh = -TT — (5) 

2n firupc Nmv^i^ 

Finally, the thirteenth column contains the ratio be- 
tween the mass outflow rate in the supershell given by 
M= 6.0X (rsh/1 kpc)(iVHi/102° cm'2)(t;sh/200 km s"^) Mq 

This implies that we convert a rest-frame UV luminosity den- 
sity into a Lya lumino sity using only the intrinsic EW given by 
IVerhamme et aL The actual star formation rates that are 

associated with these UV luminosity densities are irrelevant. 



yr~^ jVerhamme et al"]|2008h . and the UV-derived star for- 
mation rate (in Mq yr~^). 

Throughout our calcul ations of /scat we have use d the 
b-parameter simi l arly t o j Schaerer fc Verhammd (|2008l ) and 
IVerhamme et al.l (|2008l ) instead of Vth, i.e. we substituted 
Avn = hva/c in equation (Q. Below, we describe individual 
cases in more detail. 



2.3.1 MS1512-CB58 

iPett ini et al.l l|2002l ) have found that the lensed Lyman 
Break galaxy (LBG) MS1512-cB58 is surrounded by an out- 
flowing shell of gas with a column density A'hi = 7.5 x 10^" 
cm~'^, w hich is expanding at a speed of Wsh = 200 km s~^ 
(also see ISchaerer fc Verhamindl2008l ). The star formation 
rate inside the LBG, ~ 4OM0 yr~^, translates to an in- 
trinsic Lya luminosity of ~ 4-8 x 10*^ erg s~^ (with the 
uncertainty due to the unknown metallicity of the gas, see 
Schaerer 2003). We adopt the central value of Lq, = 6 x 10'*^ 
erg s~^. 

The blue dotted line in Figure [2] depicts the time evolu- 
tion of a shell with msh = 3 x IO^A^q and the above parame- 
ters of MS1512-cB58. The plot shows that the observed shell 
velocity is reached at t ~ 1.9 Myr, when the shell reaches a 
radius r = 0.21 kpc and has an observed column density of 
A^'hi ~ 7.5 X 10^" cm~^. More massive shells would reach the 
observed speed at later times, larger radii, and at a lower 
shell column density. The plot also shows that Lya pressure 
accelerates the shell to ~ 360 km s~^ after 10 Myr, at which 
point the column density has declined to ~ 4 x lO'^* cm~'^. 

Lastly, we point that the HI shell in MS1512-cB58 is 
known to contain dust, with an esti mated extinction of 
E(B - y) ~ 0.3 ([Schaerer fc Verhamme 2008, and refer- 
ences therein), which implies an ap proximate optical dept h 
at 1216 A of TD ~ 1.8 - 3.3 (e.g. (verhamme et all 120081 ). 
We argue in § [3] that despite the presence of this dust, the 
Lya radiation pressure can still be substantial (and possibly 
dominant over continuum radiation pressure). 

2.3.2 Lya Emitters m the FORS Deep Field 

IVerhamme et all (|2008l ) reproduced the observed Lya line 
profiles of 11 LBGs f rom the FORS De ep Field at 
2.8 <.z ^5 observed by iTapken et all (|2003) based on a 
simple model in which the Lya photons emitted by the 
LBGs backscatter off a single spherical outflowing shell with 
2 X lO^^cm"^ ^iVm ^7 x 10^" cm'^ 

Table [2] contains two entries for FDF 1267 because two 
different models were found to reproduce the observed line 
profile. We have not included galaxies FDF 4691 and FDF 
7539 in the table because these galaxies contained almost 
static supershells, which cannot be reproduced by our ap- 
proach, possibly because the shells in these galaxies were 
slowed down by gravity. 

For several galaxies our model can be ruled out: our 
predicted Tsh for FDF 4691 exceeds the size of the largest 
observed Lya emitting s tructure in our Un iverse, which is 
~ 150 kpc in diameter jSteidel et al.|[2OO0l ). Furthermore, 
we find ages of tsh ~ 0.1 Myr for two galaxies (FDF 5215 
and FDF 1267). While this is not physically impossible, it 
appears unlikely that we have caught star forming regions 
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Table 2. Predicted physical sizes of expanding HI supershells under the assumption that their expansion was driven predominantly by 
Lya radiation pressure. 



Galaxy A^hI fsh La r^h da m^i-^ t^h E^h Vtcrm fosc f 

name (lO^O/cm^) (km/s) (IQ-*^ £^6) (kpc) (") (IO'^Mq) (Myr) (lO^* erg) (km/s) (km/s) 



cB-58 


7.5 


200. 


6. 


0.22 


0.03 


3.3 


1.9 


1.3 


370 


134 


1.0 


0.24 




FORS Deep Field Galaxies observed by Tapken 


ct al (2007). Parameters taken from Verhamme et al (2008). 






1267 


0.2 


50. 


0.3 


5.5 


0.69 


58 


200 


1.4 


110. 


544 


0.85 


1.4 


1267 


3 


300. 


1. 


0.03 


0.004 


0.03 


0.14 


0.02 


544 


170 


1.0 


0.7 


1337 


5 


200. 


2.5 


0.12 


0.016 


0.7 


1.05 


0.27 


437 


170 


0.89 


0.1 


2384 


0.3 


150. 


5. 


7.4 


1.0 


157 


89 


35 


164 


170 


0.92 


0.4 


3389 


0.2 


150. 


1.2 


2.5 


0.37 


12 


29.5 


2.7 


210 


263 


0.87 


0.2 


4454 


0.2 


150. 


0.28 


0.63 


0.8 


0.08 


7.4 


0.2 


211 


256 


0.92 


0.3 


5215 


7 


400. 


4.4 


0.03 


0.004 


0.06 


0.1 


0.1 


780 


170 


0.65 


0.1 


5550 


5 


200. 


4.2 


0.21 


0.03 


2 


1.9 


0.84 


430 


170 


0.93 


0.1 


5812 


0.2 


150. 


2.0 


4.3 


0.67 


35 


50 


8 


196 


165 


0.90 


0.7 


6557 


0.4 


150. 


1.4 


1.7 


0.26 


11 


22 


2.5 


223. 


166 


1.0 


0.2 


4691 


0.8 


10 


18 


230 


30 


4 X 10^ 


4 X 10" 


400 




1936 


0.83 


3.9 


7539 


5.0 


25 


42 


15 


2.0 


10* 


10^ 


67 




170 


1.0 


1.9 



that early in their evolution. Indeed, these predicted ages 
fall well below the plausible age range that was derived for 
local supershells. Lastly, we predict a radius rsh = 15 kpc 
for FDF 7531. Our model therefore predicts that existing 
observations should have resolved this galaxy as a spatially 
extended Lya source (which is not observed). 

For most other galaxies, rgh ~ 0.1-10 kpc, tsh ~ 1- 
200 Myr, and E^h ~ ipS^ xQ^s ergs. For comparison, 
iMcClure-Grifflths et all (|2002h find r^h ~ 0.07-0.7 kpc, 
tsh ~ 0.9-20Myr, and E^h = 0.03-5 x 10^^ ergs for 19 galac- 
tic supershells. Our largest shells are larger than those of 
supershells that were observed in our galaxy, but are com- 
parable in size of the inferred sizes of HI outflows around 
local star burst galaxies (Mass-Hesse et al. 2003). Thus, 
there is considerable overlap in the physical properties of 
observed (extra)galactic supershells and those in our model, 
although our model does contain a few supershells that are 
significantly larger and more energetic than observed galac- 
tic supershells. This may be because we ignore gravity. In- 
deed, for these energetic shells {Enh <; a few lO^** erg), the 
maximum shell speed does not exceed the escape velocity of 
the host dark matter halo, and so gravity cannot be ignored. 

We caution that the escape velocity of the dark matter 
halo that we computed are approximate. Our escape velocity 
was derived using the relation v^sc = V^Vdrc \/ln[rvir/rsh,i], 
which assumes that the shell climbs out of the gravitational 
potential of an isothermal spheral. Although our calculated 
Vesc depends only weakly on rsh.i, we caution that devia- 
tions from our assumed simple model may change the results 
somewhat. 

Interestingly, the fudge factor / is within the range 
0.65-1.0 in all models. The small scatter in / arises because 
all shells initially have much higher HI column densities than 



The isothermal sphere model for the inner mass distribution 
of galaxies is supported observat ionally by gravitational lensing 
studies of elliptical galaxies (e.g. IWinn et al.il2003t iBolton et all 
[2OO8) and by the observed flatness of the rotation curves of spiral 
galaxies ('e.g. lBegernanlll989l :[ Sanders fc McGaughll2002h 



their observed values. Therefore, for all shells the parameter 
/scat = 1 during the early stages of their evolution. Depend- 
ing on the precise time evolution of column density and shell 
speed, /scat eventually drops below unity (for some galaxies 
this has not happened yet and /scat ~ 1, which in turn im- 
plies / = 1). Additional scatter may arise from the fact that 
some shells have predicted size that are close to the initial 
shell size assumed for our models. 

The small scatter in / implies that one can predict phys- 
ical properties of galactic supershells without numerically 
integrating equation (|3]). Instead, one may simply adopt the 
central value / ~ 0.9 and apply Eq. (O to quantities such 
as La, Nm and Vsh to predict rsh- We also find that the 
mass outflow rates in supershells reach ~ 10-100% of the 
star formation rates in their host galaxies, consistently with 
the total outflow rates that one expects theoretically (e.g. 
lErbll2008l ). 



3 DISCUSSION & CONCLUSIONS 

It is well known that scattering of Lya photons by neu- 
tral hydrogen in a single outflowing supershell around star- 
forming galaxies can naturally explain two observed phe- 
nomena: (i) the common shift of the Lya emission line to- 
wards the red relative to other nebular recombination and 
metal absorption lines, and (li) the asymmetry of the Lya 
line with emission extending well into its red wing. In this 
paper we have computed the radiation pressure that is ex- 
erted by this scattered Lya radiation on the outflowing shell. 

We have shown that for reasonable shell parameters the 
shell can be accelerated to a few hundred km s~^ after ;^10 
Myr. The Lya acceleration mechanism becomes increasingly 
inefficient as the shell speed increases (Fig. [l|, because the 
Doppler boost of the gas extends increasingly farther into 
the wing of the Lya absorption profile allowing a decreas- 
ing fraction of the Lya photons to scatter on the outflowing 
shell. Furthermore, for a shell mass that is constant in time, 
A^'hi oc , and we find that for any choice of model pa- 
rameters the shell achieves a terminal velocity that depends 



mainljifl on its total mass and the Lya luminosity of the 
central source (Fig. [Sjl . 

We have computed the physical properties of expanding 
supershells that are likely to be present in specific observed 
high-redshift {z = 2.7-5.0) galaxies, under the assumptions 
that they are driven predominantly by Lya radiation pres- 
sure and that their mass remains constant in time. We pre- 
dict supershell radii that lie in the range rsh ~ 0.1-10 kpc, 
ages in the range t^h = 1-100 Myr, and energies in the range 
-Esh = lO^'^-lO^^ ergs, in broad agreement with the proper- 
ties of local galactic supershells. We derive mass outflow 
rates in the supershells that reach ~ 10-100% of the star 
formation rates in the host galaxies, in agreement with the 
total outflow rates that one expects theoretically (e.g. lErbl 
l2008f ). 

We have found that in all models the radius of the su- 
pershell is determined uniquely by parameters such as the 
intrinsic Lya luminosity of the host galaxy (which may be 
inferred from the observed Lya line shape). La, the total 
column density of HI in the supershell, A^hi, and the shell 
velocity Vah, 

Vkpcy/ ~ V0.9/V10« A iVni A v.i, /' 

(6) 

Here / is a fudge factor that lies in the range / — 1.0-1.3. 
Our models have ignored gravity and the pressure of exerted 
by the surrounding interstellar medium. This relation holds 
up only for supershells that are driven predominantly by Lya 
radiation pressure. This relation is most accurate for those 
galaxies in which the shells are accelerated to velocities that 
exceed Vcsc- We found this to be the case in 5 out of 9 (55%, 
see Table [5J ignoring the galaxies for which our model was 
ruled out, i.e. FDF 1267b, FDF 5215, FDF4691, and FDF 
7539) of known galaxies. Hence, these shells could contribute 
to the enrichment of the intergalactic medium. 

In other cases, gravity (and external pressure) will tend 
to reduce the value of the fudge factor /. Also, if neutral 
hydrogen gas only makes up a fraction /hi of the total mass 
of the shell (see § [l]l, then our predicted value of rsh should 
be lowered by a factor of /hi. 

On the other hand, as was already mentioned in § [21 
our conservative calculations have completely ignored the 
fact that 'trapping' of Lya radiation by the optically thick 
supershell may significantly boost the Lya radiation pres- 
sure compared to what was used in this paper. Specifically, 
it was shown in Paper I that when this trapping is properly 
accounted for, the total radiation force in Eq. [T] should in- 
clude MF{vsh,Nm) in place of /scat(i^sh, A^'m). Here, Mp is 
'force-multiplication' factoi0 that depends both the column 

^ An alternative way to see why the observed quantities A''hi , "Ush i 
and La provide unique constraints on the shell mass, radius, and 
age is that these three unknown quantities relate to the observed 
quantities via three equations; 1. rrish = 47rr^jj7VHiAt»Tip and 2. 

~ tahi S'lid 3. r^ji ~ 2CTn ^sh ' found that Equations 
2. and 3. are accurate to within ~ 50%. 

^ This term derives from the (time-depen dent) force- 
multiplication function M(t) that was introduced bv lCastor et al.l 
j 19751) . as Fj-ad = M{t) ^'^^ai ^ n^^e, F^ad is the total force that 
radiation exerts on a medium, re is the total optical depth 
to electron scattering through this medium. The function 
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Figure 3. When Lyo trapping is properly accounted for, the total 
radiation force in Eq.[T] should be modified as /scat ('''sh i ^HI ) ~> 
^^('^shi ^Hl)- Here, we plot Mp (a 'force- multiplication' factor) 
as a function of the expansion velocity of the HI shell, t!sh, for 
three values of HI column density. Mp was calculated by perform- 
ing Monte-Carlo Lya radiative transfer calculation (described in 
Paper I), and under the assumption that the HI shell surrounds an 
empty cavity. Mp increases with increasing A^hi and decreasing 
Vsh- The grey line is the trajectory in the A^hi — fsh plane of the 
model of the supershell in cB-58 12. 3. 11 1. Clearly, especially at 
early times we have underestimated the radiation pressure force 
significantly, which renders our calculations conservative. 

density of HI and the velocity of the supershell. In Paper 
I, Mp was computed for a range of TVhi and v^h values, by 
performing Monte-Carlo Lya radiative transfer calculation 
for each combination of A'hi and Vsh- In these calculations, 
the HI shell surrounds an empty cavity. The results of these 
calculations are shown in Fig. [3] As evident from the plot, 
Mp greatly exceeds unity for low shell velocities and large 
HI column densities. Overplotted as the grey line is the tra- 
jectory in the TVhi — Wsh plane of the supershell in cB-58 
fij |2.3.1|l . At early times, Mp > 10, and we may have under- 
estimated the radiation pressure force significantly (as was 
mentioned in § [21 the actual force multiplication factor may 
be less than Mp, if gas and dust interior to the supershell 
prevents photons from repeatedly 'bouncing' back an forth 
between opposite sides of the HI shell) . Lya trapping could 
boost our predicted value of rsh by a factor as large as (Mp), 
where (Mp) the value of Mp, averaged over the expansion 
history of the shell. 

So far, we have assumed HI to be the only source of 
opacity in the supershell. However, supershells may also con- 
tain dust which provides an additional source of opacity; for 
example, in the models of Verhamme et al (2008), the opti- 
cal depth through dust at A = 1216 A is td = 0.0-2.0. At 
a given HI column density, we have therefore systematically 
underestimated /scat, (absorption of a Lya photon by dust 
also results in a momentum transfer of magnitude hva/c), 
which renders our calculations conservative. Furthermore, 

M(t) arises because of the contribution of numerous metal 
absorption lines to the medium's opacity, and ca n be as large 
as_Mmax(t) ~ 10'^ in the atmospheres of O-stars ll Castor et al.l 
Il975l) . 
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dusty supershells could also absorb significant amounts of 
continuum radiation, boosting the total momentum transfer 
rate even further. If the dust resides interior to the HI shell, 
then the Lya flux that impinges upon the shell is reduced 
by 6"^'''^°'. In this case, the Lya (as well as the continuum) 
radiation pressure is smaller than computed in this paper. 
However, the coupling between gas and dust makes it very 
unlikely that (Lya) radiation pressure sweeps up the gas 
surroundi ng a star forming re gion, while keeping the dust in 
place fsee iMurrav et al]|2005l . for more discussion on this). 

Of course, the existence of dust inside the HI shell re- 
duces its ability to 'trap' Lya photons, which reduces the 
value of Mf compared to that shown in Figure [3] Despite 
this reduction. Alp can still greatly exceed unity, because 
at low shell expansion velocities the majority of Lya pho- 
tons do no penetrate deeply into the HI shell upon their first 
entry. Instead, the photons mostly scatter near the surface 
of the shell, and only 'see' a fraction of the dust opacity 
(indeed, for this reason Lya radiation can escape from a 
dusty two-phased ISM, see Neufeld 1991 and Hansen & Oh 
2006). In other words, a relative small inner portion the HI 
shell can effectively trap Lya photons. We have verified this 
statement with Monte-Carlo radiative tra nsfer calculations 
that i ncluded dust at the levels inferred bv lVerhamme et al.l 
l|2008l ): we typically found Mf to be reduced by a factor 
of up to a few. In other words, even for dusty shells it is 
possible that Mf ^ 1. 

Therefore, the pressure exerted by Lya radiation alone 
can exceed the maximum possible pressure exerted by con- 
tinuum radiation (also see Paper I), which can drive the 
large gas masses (as large as a few 10^°"^*^ Mq ) out of galax- 
ies in a galactic superwind l|Murrav et aLll2005t ) . For compar- 
ison, in our model the Lya photons transfer their momentum 
onto the expanding supershell, which contains significantly 
less mass (rriah ;^1O*M0). This implies that in principle, 
both mechanisms could operate simultaneously, and that the 
neutral HI supershells may be accelerated to higher veloc- 
ities than the bulk of the ejected gas. Since in our model 
Lya radiation pressure operates only on a fraction of the 
gas, it does not provide a self-regulati ng mechanism for sta r 
formation and black hole growth as in lMurrav et al. I l|2005l ). 
However, it does provide a new way of enriching the inter- 
galactic medium with metals. 

We note that quasar s can ha ve Lya luminositiefl that 
reach La = 10^^ erg s~^ l|Fan et a l. 2006). In principle, this 
Lya luminosity could transfer a significant amount of mo- 
mentum onto neutral gas in its proximity. Because the Lya 
emission line of quasars is typically very broad (~ few thou- 
sand km s~^), /scat is significantly less than unity for the 
column densities considered in the paper (see Fig. [TJ . Fur- 
thermore, it is unclear whether gas clouds can remain neu- 
tral in close proximity to the quasar for an extended period 
of time. Even if this is the case, it remains to be determined 
whether the Lya radiation pressure can compete with the 
continuum radiation pressure. The importance of Lya radi- 
ation pressure near quasars therefore remains an open issue. 



^ Empirically, the Lya luminosity of quasars is related to their 
B-band luminosity as La ^ 0.7-Lb (Dijkstra & Wyithe, 2006). 
High-redshift quasars have been observed with Lb JjIO^^^b,© = 
4 X lO'''^ erg s'^ (see e.g Figure 1 of Dijkstra & Wyithe, 2006). 



To conclude, observations indicate that it is the kine- 
matics of HI gas surrounding star forming regions that 
mostly determines the observed propert ies of the Lya ra- 
diation (as opposed to dust content, e.g. Kunth et al.lll998l : 
lAtek et ahlboosl : iHaves et al.ll2008l : lOstlin et al.ll2008l ). This 
work suggests that -at least in some cases- the pressure ex- 
erted by the Lya photons themselves may be important in 
determining the kinematics of HI gas. This is appealing for 
two reasons: (i) the mechanism operates irrespective of the 
dust content of the HI supershells. This mechanism may 
therefore operate also in galaxies of primordial composition 
at high redshift (especially since the total Lya luminosity 
per unit star formation rate is higher, e.g. Schaerer 2003), 
and (ii) the shape and velocity offset of the observed Lya 
emission strongly suggest that momentum transfer from Lya 
photons to HI gas is actually observed to occur. 
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